The Nature of Electronic States in Atomically Thin M0S2 Field-Effect Transistors 
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We present low temperature electrical transport experiments in five field effect transistor devices 
consisting of monolayer, bilayer and trilayer M0S2 films, mechanically exfoliated onto Si/Si02 sub- 
strate. Our experiments reveal that the electronic states in all films are localized well up to the 
room temperature over the experimentally accessible range of gate voltage. This manifests in two 
dimensional (2D) variable range hopping (VRH) at high temperatures, while below ~ 30 K the 
conductivity displays oscillatory structures in gate voltage arising from resonant tunneling at the 
localized sites. From the correlation energy (To) of VRH and gate voltage dependence of conduc- 
tivity, we suggest that Coulomb potential from trapped charges in the substrate are the dominant 
source of disorder in M0S2 field effect devices, which leads to carrier localization as well. 
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MoS2-based dichalcogenides have recently been of re- 
newed interest due to the possibility of creating atomi- 
cally thin semiconductor membranes for a variety of ap- 
plications ii"— Being a layered compound, with a weak van 
der Waal interaction between the layers, M0S2 can be ex- 
foliated like graphene on insulating substrates. This has 
recently led to the fabrication of single-layer M0S2 field 
effect transistor that has a very high on-off ratio due 
to a finite bandgap.^ It has been demonstrated that the 
bandgap is indirect (w 1.2 eV) for multilayer M0S2 films 
but direct (w 1.8 eV) for single atomic layer, ■^ which may 
lead not only to low-power dissipation electronic devices, 
but also new possibilities in energy harvesting designs. 
The existence of bandgap can also have serious implica- 
tions on the charge transport and nature of disorder in 
M0S2 films, affecting its ability to screen external poten- 
tial fiuctuationSf^iS. In fact, the mobility (< 200 cm^/V.s) 
of charge carriers in single layer M0S2 devices is much 
lower than graphene , which has been attributed to the 
absence of a bandgap in pristine 2D layers of graphene. 
Our objective here is to explore the nature of disorder, 
and hence that of the electronic states, from the low tem- 
perature electrical transport in M0S2 films, when the film 
thickness is downsized from few to a single molecular 
layer. 

Often, disorder in low dimensional electron systems 
arises from extraneous sources, such as local charge dis- 
tribution that induces a random Coulomb potential on 
the electrons. This, for example, can be the remote 
dopant ions in modulation-doped III-V semiconductor,— 
or charges trapped in the substrate in case of graphene.™ 
At low carrier densities the screening of the random 
Coulomb potential becomes weak, causing carriers to lo- 
calization and/or charge distribution to become inhomo- 
geneous.S*^ This would not only have a direct impact 
on the transport, but also on the response of the system 
to various stimuli including light, stress etc. Although 



a similar charge trap-mediated transport has been sug- 
gested in thick nanoscale patches of MoS2/^ the micro- 
scopic picture is far from clear. It is also not known 
if such a picture would be valid in mono or very few 
layer M0S2 devices. In this context, our experiments 
with M0S2 devices of different thicknesses (see Table I) 
reveal that electrons are strongly localized in all cases 
which manifest in variable range hopping transport and 
an inhomogeniety in charge distribution that results in 
local transport resonances at low temperatures. We sug- 
gest that localization is probably due to strong potential 
fiuctuations induced by the randomly occupied charge 
traps that are located primarily at the MoS2-substrate 
interface. 



RESULTS AND DISCUSSION 

Devices were prepared by standard mechanical exfo- 
liation of bulk M0S2 on 300 nm Si02 on n++ doped 
silicon substrate using scotch tape technique, "'^^■-'^^ The 
fiakes were identified using optical microscope and char- 
acterized via Raman spectroscopy and atomic force mi- 
croscopy (AFM). We present detailed experiments on five 
devices (See Table I) of different film thicknesses. In 
Fig. Ic we show Raman spectra for bulk, tri and single 
layer M0S2 films. We focus on the Ej^ and Aig modes 
which have been shown to be sensitive to the number of 
atomic layers^S, The position of both modes agree well 
with recent investigation of Raman spectroscopy in thin 
exfoliated M0S2 films. The separation of Ejg and Aig 
peaks were found to be 23 cm~^, 21 cm~^ and 16 - 
18 cm~^ for tri, bi and single layer respectively. For 
further confirmation of Raman data, we determined the 
thickness of flakes using contact mode AFM. A line scan 
across the edge of a single layer fiake (Fig. Id) shows 
a step of w 0.7 nm which compares very well with the 
thickness of the single M0S2 layer (w 0.65 nm). 



TABLE I: Details of the devices: 



Device 



Number of layer 



Contact material 



Von"^ 



Device area (LxW)' 



Mobility '^ 



MoSlLa 
MoSlLb 
MoSlLc 
MoS2L 
MoS3L 



Ti/Au 


32 


4x3 


Au 


15 


2 X 2.5 


Au 


-5 


5x8 


Au 


-2 


2.8 X 2.5 


Ti/Au 


-25 


4x 16 



1 

5 
12 
20 
10 



"in V 

''both dimensions in fim 

.2, 



'^in cm /V-s near room temperature 



The electrical contacts, designed with electron beam 
lithography, consisted of thermally evaporated Ti/Au or 
Au films. The optical image of the MoS3L device is 
shown in Fig. la. All measurements were carried out 
in cryostats under high vacuum (10~^ mbar) condition. 
In all devices the gate voltage (Vbg) was applied only 
at the doped silicon backgate (see Fig. lb). Measure- 
ments were primarily two-probe current measurement us- 
ing lockin technique due to very high resistance of these 
systems, although four-probe devices were fabricated as 
well. We found, at high doping concentration, the con- 
tact resistance was negligible near room temperature, 
but increases to about half of sample resistance below 
lOOK. Detailed Ids — Vds measurements, where Ids and 
Vds are the drain-source current and bias respectively, 
were conducted to characterize the electrical contacts 
(see supplementary information). At low voltages (|Vds| 
< 300 mV), Ids — Vds at all Vbg and near room tem- 
perature, were linear for both Ti/Au and Au deposited 
samples although we have got better linear contact with 
only Au. These results bear close resemblance to the 
characteristics reported recently for high mobility M0S2 
devicesi^ As shown in the supplementary information. 
Ids — Vds characteristics become nonlinear at large Vds, 
particularly at low temperatures (T), although we at- 
tribute this to the insulating nature of the devices which 
causes the nonlinearity. The symmetric nature of Ids- 
Vds around Vds = enables us to eliminate any pos- 
sibility of Schottky contact in our operating Vds range. 
This is supported by the observed magnitude of the dif- 
ferential carrier mobility /i(= (1/C) x da/dVBa), where 
C is the gate capacitance per unit area (here 1.2 x 10"'' 
F/m2 for 300 nm SiOa), and a{= {L/W) x I /Vds) is the 
linear conductivity at low Vds- L and W are the length 
and width of the M0S2 channel. In both two and four- 
probe geometry, we obtained similar values of mobility, 
which are typical values reported for M0S2 transistors* 
on Si02 substrate (See Table I). 

In Fig. 2 we show the variation of a in MoSBL and 
MoSlLa as a function of Vbg over a range ^ 10 K to 
280 K. The conduction was achieved predominantly in 
the positive Vbg regime implying that the M0S2 films 
were intrinsically n-typeJi The doping was higher in 



MoSBL which required a negative Vbg to pinch-off com- 
pletely. Below 250 K both devices were strongly insulat- 
ing at all Vg. A weak metal-like behavior was observed at 
T > 250 K and very high doping (large positive Vbg) for 
most of the samples. Such a behavior, which was stronger 
in MoSBL, was found to be connected to a decrease in /i 
with increasing T in this regime (inset of Figs. 2a and 
2b). This can be attributed to enhanced scattering of 
electrons by phonons at high temperaturesiJ^ At low T, 
both a and ji drop rapidly with decreasing T. Such an in- 
sulating behavior was observed in multilayer nanopatches 
of M0S2 (thickness: 8 — 35 nm) as wellr^i with an appar- 
ent activated behavior of a over a rather limited range 
of T . This was explained by invoking a dense distribu- 
tion of trap states, which acted as an "impurity band" , 
although the origin or the physical location of such traps 
are unclear. 

In contrast to the nanopatcheSfii the T-dependence of 
a in our mono and trilayer M0S2 devices are not ac- 
tivated, but there are two distinct regimes in the T- 
variation in a (Figs. 3a and 3b): The high-T regime 
(T ^ 30 K), where a increases rapidly with increasing 
T, and second, the low-T regime (T ^ 30 K), where the 
variation in a weakens considerably at most Vbg in both 
devices (except for MoSlLa at low Vbc where the weak- 
ening of a sets in at higher T (Fig. 3b)). We find that in 
the high T regime, the variation of a with T can be mod- 
eled very well in terms of variable range hopping (VRH) 
transport rather than the thermally activated behaviour 
with 



a = ao(r)exp[-(To/T)^ 



(1) 



where Tq and d are correlation energy scale and dimen- 
sionality ,^^'^^ respectively, and (Tq = AT"^ with m « 
0.8 — 1. The agreement of the data to VRH transport 
with d = 2 indicates the electron transport in atomically 
thin M0S2 to occur in a wide (s> /csT) band of local- 
ized states, rather than direct excitation to conduction 
band minimum or mobility edge from the Fermi energy 
as suggested for the nanopatches. '^^ The VRH transport 
in a also results in In/i oc T^^'^ in two dimensiouj^ 
This is confirmed in the inset of Fig. Bb for the MoSlLa 
device. The magnitude of Tq decreases rapidly as Vbg^ 



or equivalently, the Fermi energy Ep, is increased. Such 
a behavior is common to strongly locahzed 2D electron 
systemSfi^ii^ and implies that Ep is located in the con- 
duction band tail. 

To understand the weakening of cr at T ^ 30 K, we 
have magnified this regime for both MoS3L (Fig. 3c) 
and MoSlLa (Fig. 3d). In both cases, the variation of 
a with Vbg becomes nonmonotonic, and displays sev- 
eral peaks which become progressively well defined as 
T is reduced. The peaks are highly reproducible, and 
stable even at T ~ 30 — 40 K, indicating that random 
fluctuations due to interference of hopping paths are un- 
likely to cause them. Resonant tunneling at the localized 
states in disordered mesoscopic semiconductors is known 
to result in strong reproducible peaks in a at low tem- 



peratures 



20^ 



In the presence of multiple overlapping 



resonances, T dependence of a weakens as observed in 
our datai^l However, confirmation of this scenario can 
be obtained by shifting the resonance peaks using finite 
Vds- For this, we focused on a small interval of Vbg 
(8 — 35 V) near pinch-off where a number of isolated res- 
onances could be identified. In the {Vbg, Vds) plane, this 
leads to diamond-like pattern in differential conductivity 
dl/dVos (Inset: Fig 3c). The occurrence of transport 
resonances indicates a rather inhomogenous charge dis- 
tribution in M0S2 films, possibly puddles of charge near 
conduction threshold, through which charging events at 
the localized states couple to the metal contacts. 

We now turn into the key issue here that concerns 
the origin of localized states in ultrathin MoS2 films. 
This requires an understanding of the origin of disor- 
der in such systems, for which we first examine the val- 
ues of Tq. However, to compare To for different devices, 
we define a device-specific reference voltage Vqn close 
to the "pinch-off" voltage in a vs Vbg curve, so that 
the difference AVbg = Vbg — Vqn is proportional to 
Ep or number density n. In Fig. 4a, we have plotted 
the variation of Tq as a function of AVbg for all the 
devices. The striking feature here is the close agree- 
ment of To in both absolute magnitude and energy over 
nearly three decades, irrespective of independent prepa- 
ration of devices, varying layer number, mobility and 
device geometry etc. This indicates a very similar dis- 
order landscape in all devices that reflects comparable 
magnitude and energy dependence of localization length 
{£,) and density of states D{E). Disorder arising from 
defects in bulk of the M0S2 films are unlikely to ex- 
plain the insensitivity of Tq to number of layers since 
screening of impurities and density of defect in bulk are 
expected to strongly influence the density of localized 
states. In stead, our data indicates a common external 
origin of disorder, such as the trapped charges in the 
substrate. This is also supported by recent transport ex- 
periments,— where higher mobility of thin M0S2 flakes 
could be achieved by changing electrostatic environment 
alone. Indeed a charge trap-induced disorder can readly 



explains the observed magnitude of Tq. To illustrate this 
we take £, as the typical size of the puddles, which for 
MoS3L can be roughly estimated to be ^ ~ 8 nm from 
the charging energy (^ 90 meV) at Vbg ~ 23 V (corre- 
sponding to AVbg ~ 48 V) (See inset of Fig. 3c). Tak- 
ing D{E) ^ 4 X 10^^ eV^^cm^^ as the typical surface 
density of charge traps at Si02 interface, ^^'^^ and using 
r, we find To ~ 6.2 x 10^ K, which is 



13.8 



^0 - kB£:^[E)D{E)' 

in good agreement to the observed magnitude from VRH 
data (Fig. 4a). 

This leads us to suggest that the physical origin of 
the localized states in ultra-thin M0S2 films is connected 
to the random potential fiuctuations from the trapped 
charges at the MoS2-Si02 interface. (See the schematic of 
Fig. 4b.) The screening of these trapped charges will be 
poor due to the large bandgap of M0S2 (unlike graphene) 
and hence can lead to a considerably long band tail. It 
is likely that the interfacial traps are randomly occu- 
pied during processing of the devices, predominantly via 
transfer of electrons from the exfoliated prestine M0S2 
layers, and subsequently form the frozen disorder land- 
scape since most experiments are conducted at low T. 

Finally to confirm the charge impurity induced disor- 
der, we have examined the nature of scattering of car- 
riers by defects at high Vbg and T so that the electron 
wavefunctions are nearly extended. If the main source 
of disorder arises from the randomly occupied interfacial 
traps, one would expect the scattering to be dominated 
by charge impurity scattering, which for two-dimensional 
electron systems with parabolic energy bands will lead 
to^ 



(T oc n^, bare Coulomb impurity (2) 

oc n, screened Coulomb impurity (3) 

In Fig. 5a, b and c, we have shown the dependence of a 
on AVbg (oc n) near room temperature, for the single, bi 
and tri layer M0S2 devices, respectively. In all monolayer 
M0S2 devices, as well as the bilayer (MoS2L) case, we find 
a oc AVgQ indicating scattering from nearly unscreened 
charged impurities. In the trilayer device (MoS3L), the 
variation in cr ^ AV^q is somewhat slower, indicating 
partial screening of the charge impurities. Assuming the 
electronic density of states to be approximately one-tenth 
of the free electron density of states at maximal dop- 
ing ('^ 5 X lO^^/cTO^) used in our experiment the Debye 
screening length in our devices can be estimated to be 
^ 1.5 — 2 nm, which is nearly three molecular layers of 
M0S2. This readily explains the bare charge impurity 
scattering in single and bi-layer M0S2, while charge im- 
purities are partially screened for the trilayer device. 

It is then natural to draw an analogy of our findings 
to other heavily researched exfoliated atomic scale tran- 
sistors in particular, graphene and topological insulators. 
The ubiquity of surface trap states probably constitutes 



a generic source of disorder in such ultra-thin field effect 
devices. Reducing substrate traps, for example, by using 
crystalline substrates such as graphene on boron nitride, 
may improve the quality of these systems considerably. 
A suspended device, as in case of graphene, could also 
lead to extremely high mobilities. 



CONCLUSION 

We have studied low temperature electrical transport 
in mono, bi and trilayer M0S2 transistor exfoliated onto 
Si/Si02 substrates. We find that the electrons in all cases 
are localized well upto room temperature at most gate 
voltages, and displays variable range hopping transport 
as temperature is lowered. We showed that the disorder 
is likely to arise from Coulomb potential of randomly dis- 
tributed charges at the MoS2-Si02 interface, and hence 
highly improved devices should be possible with appro- 
priate substrate engineering. 



METHODS 

Device Fabrication 

M0S2 flakes were exfoliated from bulk M0S2 (SPI Sup- 
plies) using scotch tape on Si02 (300 nm)/n++ Si wafer. 
To keep the disorder level comparable, the wafers were 
thoroughly cleaned by standard RCA cleaning followed 
by acetone and isopropyl alcohol cleaning in ultrasonic 
bath. The flakes with typical linear dimensions ranging 
from 2 /im to 20 /im were identified by Olympus BX51 
optical microscope. Raman spectrum are recorded us- 
ing WITEC confocal (XlOO objective) spectrometer with 
600 lines/mm grating, 514.5 nm excitation at a very low 
laser power level (less than 1 mW) to avoid any heating 
effect. The AFM measurements were carried out in con- 
tact mode with a NT-MDT NTEGRA AFM instrument. 
Ti(10 nm)/Au(40 nm) or Au(40 nm) contacts were de- 
fined using standard electron beam lithography followed 
by thermal evaporation and lift off in hot acetone. No 
Ar/H2 annealing was done in any of our devices. 
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FIG. 1: (a) Optical micrograph of a typical M0S2 device, (b) Schematic of a single layer M0S2 field effect transistor, (c) Raman 
spectrum of the bulk, tri, bi and single layer M0S2 films on Si/Si02 substrate, (d) Thickness scan along the white line across 
the boundary of the single layer M0S2 in the inset. Inset: High resolution atomic force microscopy (AFM) image of single layer 
M0S2 film on Si02 substrate. 
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FIG. 2: Conductivity a as a function of back gate voltage (Vbg) at various temperatures for (a) MoS3L with Vos = 4 mV and 
(b) MoSlLa with Yds = 100 mV. Insets show corresponding field effect mobility \i, vs temperature T at different gate voltages 
Vbg, extracted from the linear fit of a small region around a particular Vbg in the a vs Vbg graph. 
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FIG. 3: Temperature dependence of conductivity (ct) and variable range hopping (VRH) at different backgate voltages, for 
(a) MoS3L (Vds = 4 mV) and (b) MoSlLa iyos ~ 100 mV). The solid black lines are the linear fit to the data indicating 
VRH behavior in 2D M0S2 film. Inset in Fig. 3b shows variation of mobility (fi) with T~^'^ for the single layer device, (c)-(d) 
Reproducible conductance oscillations with back gate voltage (Vbg) at low temperature are shown for MoS3L and MoslLa 
respectively. Inset of Fig. 3c shows 2D map of the differential conductance dl /dVos of MoS3L as a function of back gate 
voltage {Vbg) and source-drain bias voltage (Vds) obtained at 12K. 
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FIG. 4: (a) To, extracted from VRH slope for five different devices are plotted as function of AVbg (= Vbg-Voat)- (b) 
Schematic representation of the fluctuations in conduction band of M0S2 thin films, arising due to proximity of the trapped 
charges at Si02/MoS2 interface (left) leading to the band tail and localized states (right). 
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FIG. 5: Variation of conductivity o with AVsg for (a) Single layer at 240K (diamond) and 300K (triangle, circle) (b) Bilayer 
at 300K and (c) Trilayer at 280K devices. 



